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ABSTRACT 



In airborne Anti-Submarine Warfare operations there is 
a critical requirement for maintaining an accurate relative 
plot of the sonobuoys with respect to the aircraft. This 
study proposed a method for locating sonobuoys in a pat- 
tern using ai rcraf t-to-buoy slant range information. The 
method did not use t r i angu 1 a t i on procedures and attempted 
to minimize the restrictions placed on the aircraft. The 
study showed the feasibility of the proposed methodology and 
the appro x _ i mate errors to be encountered. 
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I . I NTRODUCT l ON 



The fixed wing aircraft plays an important role in air- 
borne antisubmarine warfare (ASW) today. The primary tacti- 
cal employment of the aircraft is through the use of active 
and passive electronic sensors, the most useful of which is 
the sonobuoy. Sonobuoys are placed by the aircraft in vari- 
ous geometric patterns in the water as the tactical situa- 
tion predicts. The primary phases in the prosecution of a 
contact are search, detection, classification, localization 
and attack. Except in the rare case of a visual sighting of 
a submarine, the first four phases are almost exclusively 
conducted by the use of passive and active sonobuoys which 
transmit data, analyzed by sensor equipment on board the 
aircraft. This information is then synthesized by the air- 
craft tactical commander to yield an assessment of the tac- 
tical situation. 

One major problem in airborne ASW is the lack of accu- 
rate a i r c r a f t/sonobuoy referencing. Surface wind, sea cur- 
rents and errors in navigation quickly disperse the actual 
buoy pattern location from its assumed location. Moreover 
the buoys within the pattern frequently drift in a random 
manner, particularly when the distance between buoys is 
great. It is necessary to periodically update the actual 
buoy positions by a visual or electronic "mark on top." 

The flexibility of the aircraft is obviously reduced by this 
requirement.- In the fine localization phase where a combi- 
nation of two or more buoys in close proximity to the target 
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is providing the aircraft with fix data, it is necessary 
to fly outbound to the fix after passing over the closest 
buoy providing data. This tactic is necessary to provide 
the aircraft with final, accurate buoy location information 
just prior to the attack. Frequent missed kills are attri- 
buted to attacking a fix position generated by buoys whose 
accurate locations with respect to the aircraft were not 
known or recently updated. 

A better method for determining the location of sono- 
buoys is required. A method for locating buoy coordinates 
relative to the aircraft, remote from the buoy, would elimi- 
nate the need for "on top" repositioning. Such a method 
would also offer the additional aircraft flexibility re- 
quired for the eventual increase in the tactical load placed 
on the crew. 

A system under development [Ref. 1,2] will provide for 
a slant range measurement to the buoy. This slant range 
would be more accurate than the aircraft navigation system 
and could be converted to horizontal range knowing the air- 
craft altitude. The range to a buoy can be taken at any 
time after it has been placed in the water up to the life 
of the buoy and in any tactical environment providing the 
aircraft is at sufficient altitude to receive the signal. 

The use of the ranges provided by the buoys under the pro- 
posed system [Ref. 1] serves as the subject of this study. 

While t r i angu 1 a t i on methods obviously provide a solu- 
tion to the buoy location problem [Fig. 1] t h aircraft 



Buoy 




A/C Track 



R. = horizontal range from A/ C to buoy 

A = A/C location on track when first range 
measured 

B = A/C location on track when second range 
measured 

0 = separation angle between range measurements 



FIGURE 1 

STANDARD TR 1 ANGULAT I ON PROCEDURE 
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track, line AB, between the buoy rangings required to com- 
plete t r i a n g u 1 a t i on should be approximately equal to the 
range of the buoy. The buoy fix position is determined by 

the intersection of the range arcs to the buoy from the 

aircraft position. The errors realized in the buoy fix 
position are, for the most part, due to inaccuracies in the 
ranging system and inaccuracies in navigation from points 
A to B. If the leg AB [Fig. 1] representing the aircraft 
track is allowed to be in the order of length equal to the 
range to the buoy such that the range arcs intersect as 
close to perpendicular as possible, then the errors due to 
navigation and range measurement would be transferred to the 
buoy fix error in a nearly one-to-one ratio. If, however, 
the aircraft track, leg AB , between range measurements is 

short, inferring an angle 0 of only a few degrees, then the 

range arcs will intersect nearly parallel. The bearing 
error of the resulting buoy fix position may be extremely 
large. 

(n that the length of the aircraft track required for 
optimal error reduction in t r i angu 1 a t i on necessarily re- 
stricts the flexibility of the aircraft, this study will 
focus on a method of minimizing the error by statistical 
methods using a rapid series of rangings to the buoy over 
a short interval. It is felt that the statistical methods 
of error reduction will allow an improved solution to the 
problem over standard tri angulation. 



I I . NATURE OF THE PROBLEM 



The major portion of this study is centered around 
(l) determining a method for locating sonobuoys relative to 
the aircraft, (2) feasibility of the system with respect to 
the deployment of the aircraft in a tactical environment, 
and (3) the errors anticipated and how sensitive the pro- 
cedure is to these errors. A residual effort conducted in 
conjunction with the development of the sonobuoy location 
problem provides a conceptual solution method for deter- 
mining an accurate aircraft ground speed and true course. 

The development concepts for the ground speed and true 
course determination are not complicated and should therefore 
be adapted easily to any aircraft in the fleet today which 
would carry the associated sonobuoy ranging equipment. 

Range information fr-om the buoys will be available con- 
tinuously, sampled as desired. The slant range will be 
converted to horizontal range as described by Ellis [Ref. 1, 
Fig. A], There is no' inherent ambiguity problem such as 
mirror, images of fix positions and the solution is not 
quadrant or half plane restricted. Minor restrictions are 
placed on the aircraft for short periods of time. The air- 
craft, in most cases, will be required to fly straight and 
level at a con.tant true airspeed while sampling ranges to 
the buoys; however these intervals are qurte short, four to 
t we nty seconds as the tactical situation dictates. This 
restriction is not considered to be excessive as will be 
explained in the formal development of techniques for 
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sonobuoy location and ground speed determination. In cer- 
tain cases the measured ranges to the sonobuoys will be the 
only input required by the system, which generates a very 
appealing aspect to the solution technique. 'In other cases 
the aircraft coordinates as well as the measured range in- 
formation will be required as the only inputs. 

It should be noted here that navigational errors pro- 
ducing erroneous geographical track of the aircraft are not 
critical in the problem solution as such. The primary pro- 
blem being attacked by this study is that of locating the 
sonobuoy relative to the aircraft. If, in fact, the DR 
track is erroneous then the positions determined for the 
buoys will be erroneous; however, their positions relative 
to the aircraft will be accurate as they will be plotted 
relative to the aircraft DR position. Accurate geographical 
track and sonobuoy location information can be updated sub- 
sequent to the attack phase of the problem. Experience has 
shown that two aircraft in the process of relieving each 
other may have a relative error of up to 15 miles in geo- 
graphical DR tracks. Therefore the transfer of geographical 
positions of buoys within a pattern, is frequently unsatis- 
factory. When it is necessary to turn the problem over to 
a relieving aircraft, it will be easier to transfer a rela- 
tive picture, that is, range and bearing of all other buoys 
in the pattern from one reference buoy close to the relieving 
aircraft. 

Although the P 3 " C aircraft will far surpass- any existing 
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ASW airborne weapon system in the fleet today it is not 
anticipated that the models used currently will be phased 
out in the near future. It is important to consider devel- 
oping a system which may, with minimal modification, be 
adapted to current models of ASW aircraft. The ground speed 
and true course determination is developed with this in mind. 
Considering the contribution of locating sonobuoys from re- 
mote positions to the prosecution of a contact, it may be 
well worth while to consider the installation of a limited 
purpose computer in the current model VP aircraft vice more 
costly alternatives. 



15 



111. BUOY RANG I NG CONCEPT 



The sonobuoy has a transmitter built in for the purpose 
of transmitting to the aircraft the sensor signals received 
from the sonobuoy hydrophones. The carrier frequency asso- 
ciated with this signal is channelized in each buoy for 
buoy discrimination. A buoy ranging signal described in 
detail by Ellis [Ref. 1] will be transmitted over the exist- 
ing aircraft receivers. No new receivers or transmitters 
for the signal will be required. A modification to the 
detector section of the sonobuoy receiver will be required 
for signal separation. 

The sonobuoy ranging signal is generated by a crystal 
oscillator placed in the buoy. The method of range deter- 
mination is by comparison of the phase of the buoy signal 
with that of a standard frequency generated on board the 
aircraft. The phase difference is converted to slant range. 

The sonobuoy undergoes a warmup period in the aircraft 
prior to drop. It generates a stable frequency, f , and is 
compard with the aircraft frequency, f'. The phase dif- 
ference at this time and at zero distance is measured, A0 • 

o 

At any time, t, after the sonobuoy is dropped in the 'water 
a new comparison is made between f^ and f', and the phase 
AS^, is determined. The resulting phase difference, A0 r> 
where A 0 = A0 ^ - A0 , i s then converted to slant range; 

subsequently to horizontal range. The buoy range frequency 
has a wavelength of approximately 75 kilometers. One cycle 
difference in phase of the buoy signal from the 'reference 



is equivalent to 75 km. distance. A buoy at a range of 
10 km. and 85 km. would yield the same phase difference 
causing an ambiguity. It must be known what multiple of 
75 km. the buoy is from the aircraft. It is not anticipated 
that this ambiguity will cause erroneous ranging. 

The range errors generated by the system are covered 
in detail by Ellis [Ref. l] but one error of primary con- 
cern in the tactical employment of the system is the linear 
increase in frequency with respect to time within the buoy 
oscillator. Complete frequency stabilization is not yet a 
"state-of-the-art" realization. Testing to date with the 
oscillator currently in use in the buoy has generated an 
average drift rate comparable to a range error of 10^0 meters 
per hour. If a means was devised for the subtraction of the 
linear drift rate as a function of time, the remaining 

errors in the system could be aggregated to a classification 

2 

of normal with a mean of zero and a variance of a . This 

e 

2 

aggregation of errors being N(0,O e ) would allow a broad 

application of statistical methods for the solution of the 

error reduction problem in the tactical employment of the 

buoy range. . It is therefore assumed that in the following 

study the average linear error rate is eliminated from the 

problem and the assumption that the remaining errors are 
2 

N(0,<J £ ) will motivate the development. 
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I V . SUMMARY OF PROBLEM SOLUTION TECHNIQUES 



The following methods are developed for the solution 
of (l) So nobuoy Location, (2) Aircraft Ground Speed and ( 3 ) 
Aircraft True Course. This section provides a non-techni- 
ca 1 summary of the problem geometry, basic solution approach 
and the assumptions and restrictions governing the technique 
used in each of the three problem areas mentioned above. It 
is designed as a background preface to the more detailed 
mathematical development and should by itself acquaint the 
reader with a complete broad picture of what has been accom- 
plished. The mathematical development will be included in 
a subsequent section. 

A. S0N0BU0Y LOCATION 

It is desired to provide a system whereby the aircraft 
may, in any phase of the -tactical prosecution of a contact, 
locate any or all of the sonobuoys in the pattern with re- 
spect to the aircraft. To provide the aircraft with the 
flexibility needed in the more acute phases of the problem, 
the system must not place heavy restrictions on the maneuver- 
ability of the aircraft such as is required in standard 
t r i a n g u 1 a t i on techniques. A method whereby the aircraft 
could obtain locations from a single position and time or 
with a minimum advance leg would certainly be preferred. 

Recalling the discussion in the 1 a s t ~p a r a g r a p h of the 
introduction, an extended aircraft track, line AB [Fig. 1], 
would have to be flown in order that 0 be sufficiently large 
to reduce fix error of the buoy position. In a tactical 
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situation it may not be feasible for the aircraft to fly 
the leg AB . If the aircraft desired to maneuver between 
poirtts A and B it must do so at the possible cost of addi- 
tional navigational errors as aircraft navigation becomes 
degraded in many cases as a result of maneuvering. 

The geometry in Figure 2 was considered as a starting 
point for the development for sonobuoy location. This basic 
geometry will prevail in the description of the problem. 
Consider an aircraft with a buoy in the relative position 
shown. The aircraft is required to hold any true heading 
and airspeed for as long as it is required to measure 
ranges Rj, R 2 and R 3 [Fig. 3]. The time interval between 
each range measurement is short, approximately five seconds. 
At 200 kts. an aircraft travels 111 yds. /second. As many 
ranges as desired may be measured but more time is required 
and the additional accuracy obtained by more ranging is not 
significant as a trade off in a critical tactical situation. 
This will be a tactical decision. The study uses three 
range measurements throughout the descriptive phase. Addi- 
tional ranges were used in testing. It is assumed that in 
the short period of time the aircraft is required to fly 
straight and at a constant air speed that the coordinates 
(A.,B.) [Fig. 3 ] are accurate relative to each other. Any 
inaccuracies occuring will be absorbed in final fix error. 
The ranges R. and the aircraft coordinates (A.,B.) are in- 
puts and are assumed to be known. The buoy coordinates 
(X,Y) are unknown. 
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R. = horizontal range from A/C to buoy 

CPA = closest point of approach to buoy 
with respect to the A/C 

Q = relative bearing of buoy with respect 
to the A/C 



FIGURE 2 

BASIC TWO DIMENSIONAL PROBLEM GEOMETRY 
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(A, , B 1 ) 



(a 2 ,b 2 ) (a 3 ,b 3 ) 



Heading 



R . 

i 


= horizontal range to 


buoy as measured by 


A/C 






i = 1 ,2,- 


. . , n 


(A, ,B.) 


= A/C coordinates at 


time of measurement 


of 


range R. to buoy 


i = 1 , 2 ,. 


. . , n 


(X , Y ) 


= coordinated of buoy 


(unknown) 





FIGURE 3 

SONOBUOY LOCATION PROBLEM GEOMETRY 
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The model used to determine the (X,Y) position of the 
buoy is the method of Least Squares. The basic assumption 
fori the use of this model is that the errors are N ( 0 , a ) . 

The approach is to minimize the sum of the squares of the 
errors between the measured ranges to the buoy and the 
actual ranges. The equation for the Sum of Squares; 

n 2 

SS = £ ( 0 . — E . ) where; O.=observed value 

i = 1 

E.=expected value 
n „ 0. -E.=e. (error term) 

*- I (e .) 2 ‘ ' ' 

i = 1 

is adapted to the sonobuoy ranging problem; 

n ^ 2 

SS = £ (R.-RT) where; R.=measured range 

i = 1 

R.=expected value 

=[(X-A.) 2 H(Y-B.) 2 ] i 

r . -r r=e . 

I I I 

The approach of minimizing the sum of the error squared 
terms ' i s that of evaluating the first derivative of SS at 
zero with respect to the X and Y values desired; 

SS = l (e ) 2 = l { R - [ (X-A . ) 2 + (Y-B )V) 2 
i=i 1 i=i ' ' 1 

3 S S _ 3 S S 

3 X 3 Y °* 

This development produces a set of normal equations from 
which the values of X and Y are determined by an iterative 
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process . 



The iterative process requires that an estimate of the 
buoy angle relative to the aircraft, 0 , be supplied. This 
initial angle, 0o, provides a starting point for the pro- 
cess which by making alternately successive approximations 
of X, Y and 0 the process converges on the values of X and 
Y with the minimum error possible for the geometry. Values 
of 0o up to and including ± 15° in error of the actual 
angle 0 were tried [Fig. 4], There was little effect on 
the rapid convergence of X and Y to within acceptable values. 
In consideration of the latitude in choosing a value of 0o 
it is felt that the bearing supplied to the buoy by the DME 
indicator is satisfactory as an input value. The value of 
0o was considered to have been chosen as an estimate of the 
0 for the mit’-range measurement. 

Two studies were conducted to determine the feasibility 
of the Least Squares model for the solution of the sonobuoy 
location problem. The first study was a sensitivity analysis 
to determine the effect of the following on the convergence 
rate of the process; (l) distance to the buoy, (2) relative 
angle of the buoy, 0, (3) values of 0o as an estimate of 

and (4) number of ranges used in the solution of the pro- 
blem. 

The second study was an error analysis to examine the 
effect of range error on the final values of buoy coordi- 
nates, X and Y. Errors of varying magnitudes were applied 
to the ranges in different combinations at different relative 
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A = A/C location when horizontal range to buoy 
measured. - 

9 = actual relative angle of buoy with respect 

to A/C heading 

9 = estimated values of relative angle of buoy; 

error applied in 1° increments up to + 15° 
of actual relative angle 



FIGURE k 

RELATIONSHIP BETWEEN ESTIMATED AND ACTUAL VALUES OF RELATIVE 
BUOY ANGLE WITH RESPECT TO A/C HEADING 
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buoy angles. The results of these studies are explained in 

detail in Appendix A. 

< 

| 

B. GROUND SPEED DETERMINATION 

The development of the method of ground speed deter- 
mination is based on the concept that if an aircraft was 
approaching a buoy (0°relative bearing) which was providing 
accurate range information then the change in the horizontal 
range to the buoy with respect to time would be, exactly, 
the ground speed. This is a fast, accurate and simple way 
to determine ground speed if it could be ascertained that 
the relative bearing was exactly 0° relative. This is 
difficult to do as the DME indicator does not provide that 
accuracy. As the relative angle of the buoy increases from 
zero the first derivative of the range to the buoy with time 
becomes a function of the - cosine of the relative angle. If 
there existed a scheme for determining the change in range 
to the buoy with respect to time with the relationship of 
the relative angle then the ground speed could be determined 
for any relative angle for the buoy. 

In Figure 5 a geometry similar to that of the previous 
section is used. Consider an aircraft approaching a buoy 
as indicated. Were the aircraft to continue on this head- 
ing it would pass the buoy at its closest point of approach 
(CPA) at a point perpendicular to the aircVaft track. In 
the problem solution the aircraft chooses any buoy whose 
relative bearing is close to the aircraft heading. The air- 
craft is then required to fly straight and level at a 
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CPA = closest point of approach to buoy with 
respect to the A/C heading 

R. = horizontal ranges to the buoy measured 
at A/C locations 1, 2 and 3 

C = equal distances (assumed) generated by 
flying equal times at a constant true 
airspeed between points 1, 2 and 3 

Note: Buoy coordinates are not known, moreover 

are not required for this problem. 



FIGURE 5 

GROUND SPEED DETERMINATION PROBLEM GEOMETRY 
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constant true air speed for a time interval sufficient to 
measure from three to eight ranges to the buoy in question. 
The time between each measurement is short, two to four 
seconds. The time between each range, C, should be exactly 
equal. This could be done by electronic time gating. This 
time, C, also represents equal distance between range mea- 
surements. It is assumed that for short periods of time and 
by electronic time gating, these distances will be equal. 

Any errors are absorbed in the total ground speed error. 

The actual or estimated position of the buoy is not 
required as an input to the problem. Also it is not re- 
quired to know the heading or airspeed of the aircraft, 
only that they are constant. The evaluation of C deter- 
mines the ground speed. The use of the Pythagorean theorem 
produces three equations for the three ranges measured; 



The magnitude of the ranges only are known and not the direc- 
tion. There' are three equations and three unknowns; X, Y 
and C. The values for C are evaluated from the above equa- 
tions for the three ranges; 



By a similar development C is determined for any- number of 



2 2 2 

(X+2C) + Y = Ri for range measured at 1 

2 2 2 

( X+ C) + Y = R 2 for range measured at 2 

2 2 2 

X + Y = R 3 for range measured at 3 




C 




2 
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ranges measured. Up to eight ranges are used in the de- 
tailed mathematical development included later. 

Two studies were conducted to determine the feasibility 
of the ground speed determination method. The first study 
was a sensitivity analysis to determine the effect on ground 
speed error by varying the following; (l) number of ranges 
used, (2) time between range measurement, (3) ranges to the 
buoy and (4) the relative bearing of the buoy. 

The second study was an error analysis on the ground 
speed as affected by errors in the measured ranges. The re- 
sults of these studies are combined in Appendix B. j-j f 7 , 

This method for determining ground speed is expected to 
have limited application, however, because of the simplicity 
of concept it should be quite reliable to supplement break- 
down of more sophisticated inertial and dopplar navigational 
systems. Frequently, in fine localization, a timed run out 
from some reference to a fix position for attack is required. 
It is felt that this ground speed determination method can 
be utilized on the final leg to more accurately determine 
the ground speed and hence a more accurate time to fix on 
that leg. The aircraft will maintain a true course and air- 
speed to that fix position which meets the requirements for 
the solution. 

C. TRUE COURSE DETERMINATION 

/ 

It was felt that the use of the range measurement to a 
buoy as the only input could provide additional information, 
if applied properly, in evaluating a "drift vector". Consider 
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in Figure 6-a, an aircraft flying a true heading of due 
North. This true heading is held, at a constant true air 
speed (TAS), for a time t^. Vector represents the true 

air speed times the time flown, t^. The aircraft then 
turns to a true heading of East and holds this heading, at 
the same TAS, for a time t ^ where t ^ = t ^ . While on each 
of these headings the aircraft determines ground speed by 
methods described in the previous section. A different 
buoy is used for ground speed determination for each of the 
two headings. If a wind from the Southeast were affecting 
flight, the ground speed for the northerly leg would be 
greater than the TAS and the ground speed for the easterly 
leg would be less than the TAS, [Fig. 6-b]. The triangle 
shows the second leg geometry without the advance 
and transfer of the aircraft in the turn. The magnitude 
only and not the direction of these ground speed vectors is 
known. It should be noted that any initial true heading 
and TAS could be flown and a subsequent right or left turn 
could be executed. Maximum accuracy of the drift vector is 
obtained when the subsequent heading is perpendicular to the 
initial head i ng . 

In the development for the determination of a drift 
vector, the advance and transfer of the aircraft was elimi- 
nated in further diagrams and does not alter the geometric 
accuracy of the development. The actual path of the air- 
craft, less -the advance and transfer, [Fig. 7 - a], is repre- 
sented by the vectors G.. and G,. whereas the vectors D,. and 

N E N 
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Advance and transfer 




D„ = t. 



D r = t 



W = t 



N 



W„ = t, 



G„ = t, 



G„ = t, 




where t = t r 
N E 



°N " °E 



W 1 " W 2 



(Windspeed + additional errors) 

(Windspeed + additional errors) 
(Groundspeed for North leg) magnitude only 
(Groundspeed for East leg) magnitude only 



FIGURE 6 

TRUE COURSE DETERMINATION PROBLEM' GEOMETRY 
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W = drift vector = wind vector plus additional error 

components 



Known Values: 

A/C assumed position coordinates, (a,b), (c,d) arrd (e,f) 

vectors and D^-, both magnitude and direction 

vectors and G^, magnitude only 

Unknown Values: 

vectors G^ and Gr, direction only 
coordinates (X , Y } 

drift vector, W, both magnitude and dt-rection 



FIGURE 7 

TRUE COURSE DETERMINATION SOLUTION GEOMETRY 
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represent the "no wind" track. The geometry of the pro- 
blem was altered for the purpose of the mathematical solu- 
tion technique. The triangle was shifted [Fig. 7"b] 

such that coincided with since . The vector 

did not change in orientation thus the relationship be- 
tween the two triangles, and G^W^D^, was not changed. 

The resulting form revealed that the vectors G^ and G^ had 
the same terminal point. The known and unknown values in 
the figure are listed at the bottom of Figure 7 < 

A pair of quadratic equations were used to solve for 
the terminal point, (X,Y); 

(X-a ) 2 + CY-b) 2 = G 2 
(X-c) 2 + (Y-d) 2 = G 2 

The magnitude of the drift vector, W, was then solved > as 
with the determination of the coordinates (X, Y) both the 
coordinates of its origin and terminal point were known. 

W = (X-e) 2 + (Y-f ) 2 

The resulting drift vector was then used in conjunction with 
TAS and true' heading to determine an actual course and 
ground speed . 
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V . SONOBUOY LOCATION TECHNIQUES 



From the description of the basic sum of squares model 
in section IV, the sum of squares equation for the parti- 
cular geometry in Figure 3 was obtained; 

SS = I l {R. - [(X-A. ) 2 +(Y-B. ) V} 2 
i = 1 

In further development let 



x. = X-A . 

i i 



y. = Y-B. 

i i 



(1) 

( 2 ) 



then 



1 



2 2 i 2 



SS = - > [R. -(x.+y . )*] 

n . L , i i ' i 

i = 1 



(3) 



Equation (3) was then differentiated with respect to X and 

v . , , . 3 S S 3SS 

Y and evaluated at zero, i . e . , ^ = 0 



} iK i - ( x% Y ! , i, ( -i ) u x .Kvyjr* - o 

i = 1 



1 



- 1 l 

n L* 



R.x. 



i i 



“ 1 [ x. » 0 

2 -i - n 1=1 1 



i = 1 



(x.+Vj ) 2 






Let (a) be 



1 l 

n C 



R.x. 

i i 



i =1 2 2 i 

(x. +y . ) 

i i 
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and let (b) be 



1 I * 

n , L . 



i = 1 



Each term of the expression (a) was multiplied, both numera- 
tor and denominator, by — 

X • 



1 l 

n . L 



R.x. (7 ) 
1 1 x i 



R . 



R 



: _ 1 1 2 2 

( L ) (x.+y.) 5 

X . I 1 I 

I 



- = - I — - — = 1 y 

i n i-1 y 2 , "i = i 



I 



d + -;) 



( 1 + 1 a n 2 0 . ) ^ 

i 



where 0. is the relative angle of the buoy at the time the 
range R. was measured, and 



tan 0 . = 



y. Y-B. 

' i i 



i x . X -A . 



( 5 ) 



Thus (a ) become s ; 



1 n 

1 l 

n . L . 



R. 



( Se c 2 0 . ) 1 



= — J R . cos0 

n *-» i 



i = 1 



i i 



The second term, (b) was expanded; 



< n 1 n v n . n ■« n 

1 y x .= - i (x - a . ) = y - 1 y A .= x- - 1 - y 

n i = l 1 n i = 1 ' n i = 1 n i = 1 1 n i = 1 



A 



The two terms (a) and (b) , were then combined; 
- n n 

- y R.COS0. -x + - y A. = 0 

n . L . i i n . L . i 

i=l i=1 



1 n i ■■ 

X = - y R.COS0. - - y A. 
n j £; j i i n , ~ , i 



n 

I 

i = 1 



( 6 ) 



3 ^ 



A similar development for 



ass 

av 



0 yielded 



n 



n 



Y = - 7 R.sinO. - - 7 B. 

n . L , i i n . , i 

i=1 i=1 



( 7 ) 



The unknowns X and Y could not be solved directly as 0. is 
a function of X and Y (see equation 5); therefore, an initial 
value common for all the 0. was estimated, namely 0 Q , and 
initial values for X and Y, i.e., X Q and Y Q were determined 
using 0 Q for each 0. in equations 6 and 7- These values of 
X 0 and Y q In turn were used to update the values for 0. , 
namely; 

Y - B . 

tan-. = x ° _ A | 'Y'. i = 1 ,2 , . . . ,n 

These values for the 0. were then substituted in the equa- 
tions 6 and 7 to determine X x and Yj. This process was 
repeated until the convergence of X and Y to within accept- 
able values was attained. 

The foregoing iterative process proved to be lengthy. 

The angles 0. were quite insensitive to small changes in 
the positions of X and Y through the iterative process. It 
was felt that an application of Newton's method in two vari- 
ables might provide a more rapid convergence. The method 
involved the expansion of the functions (as determined by 
equations 6 and 7 using 0 Q ) in Taylor's Series in two 
variables about X 0 and Y Q . The expansion would terminate 
after the linear terms and the remainder ignored. The 
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The resulting equations would be solved for X and Y and the 
solutions used for an updated expansion. This process would 
be repeated until the values of X and Y were within accept- 
able limits. 

The sine and cosine of the relative angles 6. as ap- 
pearing in equations 6 and 7 were developed in terms of x. 
and y. (see equations 1 and 2); 



cos 0 . = 

i 



1 

(sec 6 . ) r 



( 1 +tan 0 . ) 



(l+Yj/xj)* 



x . 

i 



/ ^ 2 » 4- 

(x.+y . ) 



- 1 2-?- p o J- 

1 (esc 0.) 5 (l+cot e.) 2 (l+x!/y.) 2 

I I I I 



I 

2 2 T~ 

(x.+y.) 5 

i i 



The equations 6 and 7 were then cast into the homogeneous 
forms below; 



gi (X.Y) - X - 1 l (R ' - A.} = 0 

1=1 (x.+y.) 

i i i 

9 2 (X , Y ) - Y - I J ( R ,-.- p v . - B.} - 0 

1=1 (x.+y.) 



The form of the Taylor expansion is 



g(X,Y)=g(X 0 ,Yo)+g Y (Xo ,Y 0 ) (X-X Q )+g(X Q ,Yo) (Y-Y 0 )+Remainder 

x y 
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Specifically, expanded; 



! gi (X,Y) = Xc-I l {R 



X . 



i = 1 ' (x?+y?) 2 

I 7 I 






+ 4 } R i , - 0 (8) 



. h if 2 2 \ V2 

1=1 Cx . +y . } 

1 7 1 



and g 2 expanded; 



g 2 ( x » Y ) = Y 0 -T~ l (R, 



2 2 



i = 1 ( x . + y . ) 

l 1 I 



T “ B i }+[ ?r J, R i 



x . y ♦ 

2 ‘ 2 I 3 > -](X-X 0 ) 

i = 1 ' U*Y.) ' 2 



+ [i -- i r. 

n ^ 1 



2 
X . 



1 V r* I 

1=1 ' (x . +y . ) /2 



n .* "i / 2 2 ^ KY-Y 0 ) = 0 (9) 



with x. and y. evaluated at X 0 , Y 0 . 

Equating (8) and (9) to zero yielded the following system 
of linear equations; 

2 

1 n y • 1 n x . y . v 

t>4 l . a Y S/ . )X + [J- l R j : 2 



0 • “ 1 ' ( * 2 \ V 2 

=1 (x.+y.) 



i=1 '(x.+y!) 3/2 

1 ' 1 



Y 0 n x . y . 



= - l R. 

n 1 



I ' I 



x„ n 



n • u 1 > / 2 , 2 \ V 2 n . , R i 7 2 . — 2 \ 3 /2 

1=1 (x.+y,) z 1=1 (x.+y.) 



+ 1 l 

n 



n x . 



2 2,1 



n i = 1 Cx . + y . ) 



1 I A - 

n . L , 1 



i = 1 
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1 n x . y . 

1 V n I I 



[- I R. 

n . L , i 



. - , ./22.3/ 2 

1 = 1 ( x j + y . ) 



JX - (1-i } R, 



2 
X . 

I 



i = 1 1 ( x 2 . + y 2 . ) 3/ ^ 2 

i ' i ' 



]Y 



X n x.y. Y „ n 

_° Y r ' ' __° y 

n . Z i / 2 2 3/2 n . Z 

i = 1 (x.+ Y.) 

i 



2 
X . 

I 



t , 2 2,3/, 

1 = 1 (x.+y . /2 

i ' i 



1 I 

n 



Q y : 



+ - l B 
n i = 1 



n 



n . . / 2 , 2 \i n . , i 

i - 1 (x . +y . ) 5 • 

i ' i 



Cramer's rule was applied to solve for X and Y. The coef- 
ficient matrix, H, of X and Y, 



H = 



D ! °3 



D 2 D 4 



1=1 l R 

n L* 



- I R.- 

n L* i 



x.y. 

i 7 i 



n . 1 i / 2 ^ . 3/2 n i / 2 2 . 3/2 

1=1 (x.+y.) z i=1 (x.+y.) 



« n x.y. 

- ; 1 . 

n L* 1 



n . ^ , 1 / 2 2 3/ 2 

i=1 (x.+y.) 

1 7 1 



1 ' 1 



1 y r . — _ , , 

" i=1 '(x'+y*)* 



n 



2 
x . 

1 



and the vector of constants, 



Y n n 

- I R. 

n Z-* 1 



x : y 



x n 
0 



- I R- 

n 1 



i . 1 



n x 



i . 1 



,, ,, . — y ! h — Y 

n - i ' r 2 , 2 \ 3 /2 n . L i , 2 2.3/2 n . L . , 2 2.1 n . ^ , 

i = 1 (x.+y.) i=1 (x.+y.) i = 1(x.+y-.) 5 i = 1 



1 1 



1 1 



1 1 



x„n x.y. Y„n x 2 ,n y. ,n 

_° Y R ■ 1 __° Y R ' + I Y 1 +1 Y 

n . i , 2 2 ■> V2 n . Z . i/ 2 2 . 3/ 2 n .. / 2 2 , n ,Z 

i=1 (x.+y.) i=1 (x.+y.) i=1(x.+y.)5 1 = 



1 * 1 



1 * 1 
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were applied in the following manner; 




D2 D 4 



( 10 ) 



( 11 ) 



Equations 10 and 11 replaced equations 6 and 7 in the itera- 
tive process. This method proved highly successful as a 
rapid convergence process. The results of the analysis of 
this process are represented by the figures in Appendix A. 

With the success of the second method an attempt was 
made to interpret the distinctive improvement. The coef- 

P 

ficient matrix, H, was broken down and analyzed; 

2 _ 



- 1 R- 

n i = i ' 



(x j +y J ) 



]_ 

n 



n x . y . 

V R 1 1 

L * J 2 2 . 3 a, 

i = 1 (x.+y . ) /2 



H = 



! n x . y . 

- I R. — - 

n i=1 ' (x!+y-) 



3/ 2 



K .I/iTTrrr 

i=1 (x.+y.) 



3/2 
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The elements of the second term were then factored: 



I 



n 





R 

— br( sin 

(x.+y . ) ^ 



V 



J, 2 (sine cose,) 

1=1 (x.+y. )? 



R. 



I 



2 2 
(xT+yT)* 

i ' i 



-(sin0.cos0.) -j 

i i i L 



R. 



1-1 (x j +y? ) i 



(cos 



e i ) 



The matrix below extracted from the second term of 
the coefficient matrix has a determinant whose value is 
equal to zero for any values of 0.; 



2 

-sin 0. 

i 

s i n 0 . co s . 

i i 



s i n 0 . cos 0 . 

1 i 

2 



“cos 0 . 



If, in fact, there were n_o errors between the measured 

2 2 i 

ranges and the actual ranges then R. = (x. + y.)? and the 

R. ill 

factored term, , would equal one. This implies 

(x j + y j ) * 

that if there were no error between the measured and actual 
ranges, the second term in the matrix H would go to zero 
leaving; 



H 




0 



0 

1 



It appears that the existing error between the actual and 
measured ranges causes variations in the value of the second 
term in the -coefficient matrix H and therefore provides the 
driving force of the convergence process. 



VI. GROUND SPEED DETERMINATION 



From the problem description in Section IV, the mathe- 
matical model below for ground speed determination was de- 
veloped. Figure 5 provided the geometry for the problem. 

The particular geometry in the figure indicates three range 
measurements taken at equal intervals. Ground speed, using 
from three to eight ranges measurements was developed in 
this study. The distances Cjcalculated below are inputs to 
the ground speed determination equations 18 and 19 at the 
end of this-section where i = 3 , 4 , . . . , 8 represents the num- 
ber of measured ranges utilized in the calculation. 

For the three range geometry the known values are R , 

R 2 and R It is also assumed that the time between range 
measurements is the same in each case. The unknown values 
are X, Y and C. The thre_e equations from which C was deter- 
mined are applications of the Pythagorean theorem. 



(X + 2 C ) 2 + Y = R 2 



.2 .2 2 2 2 

(X + C) + Y = R X + 2 CX + 

2 2 2 
X + Y + R 3 X 



+ 4CX + 4c + Y = R (a) 

l 

2 2 2 
c + Y = R„ (b) 



2 

2 2 

+ V =8, (c) 



Equations (a) plus (c) mi nus 2(b) ; 



X 2 + 4CX + *»C 2 + Y 2 = 



X + Y = R 

- 2 X 2 - 4CX - 2 C 2 -2Y 2 = - 2 R 



2 C 



— Rj “ 2 R 2 Rg 



c 3 -■ 



2 2 2 
R j - 2 R 2 + R 3 



( 12 ) 
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* 



The development for four ranges; 



(X + 


. 2 2 

3 C ) + Y 


= 


2 

r i * 


2 

X 


+ 


6CX 


+ 


2 

9C 


+ 


Y 2 


= R ? 


(a) 


(X + 


2 C ) 2 + Y 2 




- 


X 2 


+ 


ACX 


+ 


Ac. 2 


+ 


Y 2 


" R 2 


(b) 


(X + 


C) + Y 


= 




X 2 


+ 


2CX 


+ 


c 2 


+ 


2 

Y 


2 

= R 3 


(c) 


X 2 


CM 

+ 


= 




X 2 










+ 


Y 2 


- R 4 


(d) 


s (a) 


minus (b ) 


minus 


(c) 


plus 


Cd) 


> 










2 

X + 


6CX + 9C 2 


+ 


II 

CM 


r : 



















2 2 2 2 

-X - ACX - AC - Y = R 2 



-X - 2 CX - 

..2 



2 2 2 
C - Y = R 3 

,2 _ 2 
+ Y = Ri, 



AC‘ 



Cu = 



= Rj - R 2 -R3 + Rj 

, “2 ~2 “2 ~ 2 ~ 
f R 1 ‘ R 2 “ R 3 + R 4 



(13) 



The development for five through eight ranges was omitted 
but was similar to the two cases above; 



C = 
5 



, 2 22 

|2Rj - 3R 2 + R, 



2 2 
R 4 + R 5 



OA) 



c = 

6 



i 2 2 2 2 2 2 

2R, - 2 R - R + Rf " R. + R: 

1 2 3 4 5 6 



12 



(15) 



C 7 = 



2 2 2 2 2 2 
R + R - 2R, + R c - R + R„ 

23 4 5 6 7 



1 2 



06 ) 



C 8 - 



R ^ + R ! " R f " R c + R 2 - R 2 + R 2 

2 3 4 5 6 7 8 

~Te 



(17) 



The value C. is a distance in units of nautical miles 

1 



or yards in each of the cases above. 



A2 



The conversion to ground speed was easily accomplished; (l) 
for C. in nautical miles and (2) for C. in yards; 



(1) GS 



C i 3600 (sec) 
T I ME ( sec ) 1 (hr) 



3600 • C 
TIME 



■ l knots 



( 18 ) 



(2) GS 



^ i S ^ 3600 (sec) 1 ( nm) 

T I ME (sec) 1 ( h r ) 2000 (yds) 



1 8 • Ct 
T I ME 



knots 



09 ) 



where TIME is the time interval in seconds between each 



range measurement. 



VII. TRUE COURSE DETERMINATION 



The true course determination used the ground speed 
calculated in the previous section. Figures 6 and 7 of 
t i on V provided the geometry for the problem and Figure 
is annotated for the following mathematical development. 
The coordinates of X and Y were first determined. This 
done by the simultaneous solution of a pair of quadratic 
equations. 

(X-a) 2 + (V-b ) 2 - G 2 
(X-c) 2 + (Y-d) 2 - G 2 



Expanding these equations; 



2 2 2 
X - 2 a X + a + Y - 


2 b Y + b 2 


= 


g n «•) 




2 2 2 
X - 2cX + C + Y - 


2d Y + d 2 


= 


G 2 (b) 




Subtract (b) from (a); 










2 2 

2(c-a)X + a _ c + 


2 (d - b ) Y + 


b 2 


2 2 

- d = g n - 


g e 


Solve for Y in terms of X 


9 








2 (d -b) Y = 2 (a -c ) X + 


g n - g e 


H" 


c 2 + d 2 - a 2 


« b 2 


y _ (a-c). y 


G N ’ G E 


4- 


c 2 + d 2 - a 2 


- b 2 


Y (d-b) X + 






2 (d-b) 





Let 

R 

Then 



a - c 

d - b 



and 




2 2 2 2 ? 
G ^ +c +d -a -b 



2 (d'b) 



kk 



a s 
Sec 
7“b 

wa s 



Y 



R ' X + S 



( 20 ) 



Substitute the value for Y in equation «(b) above; 

X 2 - 2cX + c 2 + (R-X+S) 2 - 2d (R • X+S ) + d 2 = G^ 2 

X 2 - 2cX + c 2 + R 2 X 2 + 2 RS X + S 2 - 2dRX“ - 2dS 

2 2 

+ d - G E = 0 



Col 1 ect i ng terms ; 

2 2 2 2 2 2 
(1+R )X + (2RS-2c-2dR)X + (c +d +S -2dS~G E ) 

This is a quadratic in X; 

v _ -B + B 2 AC 
X 2 A 



0 



( 21 ) 



where 

A = (1+R 2 ) 

B = (2RS-2c-2dR) 

C = (c 2 +d 2 +S 2 -2dS-G E 2 ) 

Equation (18) is a quadratic form providing two possible 
solutions. The values of the coordinates X and Y determined 
by both solutions of X provide points which are mirror images 
on each side of a line joining the points (a,b) and (c,d) 
in Figure 7 - b. The choice of coordinates (X,Y) depends on 
the initial heading of the aircraft and whether a right or 
left turn is executed on the second leg. 

From the coordinates (X,Y) the drift vector was then 
determined. The coordinates (e,f) are known. 

W = ^|( X -e ) 2 + (Y-f ) 2 - (22) 

The drift vector was then applied to the terminal heading 
to determine the true course of the aircraft for., accurate 
track determination. 
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VIII. CONCLUSIONS 



^ The study showed that the method of determining sono- 

i 

buoy location relative to the aircraft with .consecutive, 
shor t - i n terva 1 range measurement techniques is feasible. 
Moreover the application of the principle of the first 
derivative of the range with respect to time of a buoy near 
the aircraft heading as a means for determining ground speed 
is accurate and relatively insensitive to range errors, num- 
ber of ranges used and the relative angle of the buoy. 

An increase in the number of consecutive ranges mea- 
sured in excess of three, reduces the error in buoy location 
but the reduction is slight. The aircraft should not sac- 
rifice flexibility to acquire additional range measurements 
unless tactics permit. 

In the sensitivity analysis of the sonobuoy location 
technique, the convergence process appeared to break down 
for acute relative buoy angles where the estimate of the 
relative angle lay between the actual relative angle and 
the ai'rcraft heading. Before implementation of this system 
further research should be conducted in this area to pos- 
sibly provide a 360° zone of reliable information on buoy 
location. The model, at present, can provide 360° of re- 
liable information. For acute relative buoy angles, with 
respect to the aircraft heading, however, the estimated 
angle must be greater than the actual angle. 

It is felt that the utilization of buoy ranging in-, 
formation for the location of sonobuoys, apart from 



t r i angu 1 a t i on procedures, has a broad field of application, 
yet uncovered. A preliminary study was conducted to deter- 
mine relative buoy positions independent of aircraft navi- 
gation. Ranges to all buoys in the pattern were measured 

from a single aircraft position. For a pattern consisting 

N ! 

of N buoys there were ! possible two-buoy combi- 

(2!) (N -2 ) ! 

nations. The actual buoy positions were not known thus the 

assumed positions were used. The ranges were plotted in 

reverse from the assumed positions of the buoys generating 
N ! 

1 aircraft fix positions determined by the ranges 

(2!) (N-2) ! 

from each two-buoy combination. From these resulting fixes 
an average position was determined by least squares method. 
An iterative process was then developed whereby the assumed 
buoy positions were adjusted such that the reverse plot- 

M i 

ting of the ranges wou 1 d .generate 1 aircraft fixes 

(2!) (N-2) ! 

that were coincident. This did not produce a unique set 
of final assumed positions for the buoys but did establish 
a locus for each buoy in the pattern relative to the air- 
craft _ pos i t i on . It is felt that further study in this area 
would provide, possibly, a more versatile use of sonobuoy 
ranging. 



hi 



APPENDIX A, SONOBUOV LOCATION ERROR ANALYSIS 



The first study conducted with respect to sonobuoy 
location was to determine the response of the speed of 
convergence of the iterative process as affected by the var- 
iance of ; 

(1) The relative angle of the buoy with respect to the 
aircraft heading, 6. 

(2) Varying estimates, 0 q , of the relative angle. 

(3) Range from the aircraft to the buoy. 

Figures 8 through 11 show the effect on the number of itera- 
tions required to converge by variance in the estimate of 
the relative angle of the buoy. The actual relative angle 
ranged from 10° to 80° in ten degree increments. The esti- 
mate of the relative angle, in all but one case, ranged 
from -15° to +15° in error of the actual relative angle in 
increments of one degree. For the case where the relative 
buoy angle equals 10° the angle estimate, 0 q , started at 0° 
which was -10° in error of the actual angle. The convergence 
is quite rapid for all estimates, 0 , particularly when 

9q w 0. The convergence process breaks down when the buoy 
is close to the heading of the aircraft and an estimate 0 q , 
between the actual relative angle and the aircraft heading 
is chosen. This can be seen in Figure 8. It should be 
noted that for any relative angle, 0, however small, if an 
estimate, 0 q , is chosen which is less acute than Q, the 
process provides rapid convergence. Data for the effect 
on the speed of convergence for variance in buoy range is 



not included. The effect was relatively insignificant. 

Buoy ranges at 30 nm. and 5 nm. were tested with at most 
one: iteration difference tin favor of the 30 nm. range) in 
the convergence. The data in Figures 8 through 11 repre- 
sent a buoy range of approximately 15 nm. and a range error 
of zero. The convergence process stopped when the buoy 
location error was within .5% of the range to the buoy. 

The second study incorporated a range error to deter- 
mine its effect on final buoy location at the termination 
of the convergence process. Relative buoy angles of 20°, 
^5° and 80° were investigated. Figures 12, 13 and 1 ^ show 
the total buoy location error in yards versus the appli- 
cation of the range error in the following manner to the 
three ranges ; 



R R R 

0 0 0 

+ + + 



+ 



+ 



no error applied 
all errors added 
all errors subtracted 
other combinations 
other combinations 



There appears to be an approximate one-to-one relationship 
of the average buoy location error to the range error ap- 
pl ied which is appealing with respect to the small angles 
utilized between range measurements. 
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Iterations for convergence Iterations for convergence 
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FIGURE 8 

NUMBER OF ITERATIONS REQUIRED FOR CONVERGENCE VERSUS 
ESTIMATE OF RELATIVE BUOY ANGLE THETA 
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FIGURE 9 

NUMBER OF ITERATIONS REQUIRED FOR CONVERGENC-E VERSUS 
ESTIMATE OF RELATIVE BUOY ANGLE THETA 
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FIGURE 10 

NUMBER OF ITERATIONS REQUIRED FOR CONVERGENCE VERSUS 
ESTIMATE OF RELATIVE BUOY ANGLE THETA 
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Iterations for convergence Iterations for convergence 





Estimated value of relative buoy angle in degrees 

FIGURE 11 

NUMBER OF ITERATIONS REQUIRED FOR CONVERGENT VERSUS. 
ESTIMATE OF RELATIVE BUOY ANGLE THETA 
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Error in buoy location in yards 




Method of application of range error 



FIGURE 12 

ERROR IN BUOY LOCATION IN YARDS VERSUS METHOD O.F APPLICATION 

RELATIVE BUOY ANGLE 20° 






Error in buoy location in yards 




FIGURE 13 

ERROR IN BUOY LOCATION IN YARDS VERSUS METHOD OF APPLICATION 

RELATIVE BUOY ANGLE A5° 
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Error in buoy location in yards 
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FIGURE 14 

ERROR IN BUOY LOCATION IN YARDS VERSUS METHOD OF APPLICATION 

RELATIVE BUOY ANGLE 80° 
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APPENDIX B, GROUND SPEED DETERMINATION ERROR ANALYSIS 



The results of the study conducted on the ground speed 
determination are pictured in Figures 15 through 20 . They 
represent the ground speed error in knots versus the range 
error applied, under the varying conditions listed below; 

(1) Number of ranges measured (3, A,.. .,8). 

(2) Time between range measurements. 

(3) Relative bearing of buoy to aircraft heading. 

In no case does the ground speed exceed 1.75 knots. 

/ 

The case of-the relative buoy angle of zero was not tested. 
It was felt that this was a trivial case. The change in 
horizontal range with respect to time of a buoy dead ahead 
was the ground speed exactly. It was desired to examine 
the accuracy of the method of application of first deriva- 
tive of range with respecj: to time as the relative buoy 
angle increased from zero. 
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Ground speed error in knots Ground speed error in knots 




Range error in nautical miles 




Range error in nautical miles 
FIGURE 15 

GROUND SPEED ERROR IN KNOTS VERSUS ERROR IN RANGE DUE TO 
DRIFT IN BUOY FREQUENCY WITH RESPECT TO TIME - 3 RANGES 
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Ground speed error in knots Ground speed error In knots 




Range error in nautical miles 




Range error in nautical miles 
FIGURE 16 

GROUND SPEED ERROR IN KNOTS VERSUS ERROR IN RANGE DUE' TO 
DRIFT IN BUOY FREQUENCY WITH RESPECT TO TIME - b RANGES 



59 



Ground speed error in knots Ground speed error in knots 




Renge error in nautical miles 




Range error in nautical miles 
FIGURE 17 

GROUND SPEED ERROR IN KNOTS VERSUS ERROR IN RANGE DUE TO 
DRIFT IN BUOY FREQUENCY WITH RESPECT TO TIME - 5 RANGES 
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Ground speed error in knots Ground speed error in knots 




Range error in nautical miles 




Range error in nautical miles 
FIGURE 18 

GROUND SPEED ERROR IN KNOTS VERSUS ERROR IN RANGE DUE TO 
DRIFT IN BUOY FREQUENCY WITH RESPECT TO TIME - 6 RANGES 
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Ground speed error in knots Ground speed error in knots- 
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Range error in nautical miles 
FIGURE 19 

GROUND SPEED ERROR IN KNOTS VERSUS ERROR IN RANGE DUE TO 
DRIFT IN BUOY FREQUENCY WITH RESPECT TO TIME - 7 RANGES 
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Ground speed error in knots Ground speed error in knots 




Range error in nautical miles 




Range error in nautical miles 
FIGURE 20 

GROUND SPEED ERROR IN KNOTS VERSUS ERROR IN RANGE DUE' TO 
DRIFT IN BUOY FREQUENCY WITH RESPECT TO TIME - 8 RANGES 
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SONOBUOY LOCATION PROGRAM 






SONOBUOY LOCATION SIMULATION PROGRAM 



START 
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GROUND SPEED DETERMINATION PROGRAM 
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SONOBUOY LOCATION PROGRAMo W I TH . MOD I F I C AT I ONS , TO BE USEO 
FOR ACTUAL DETERMINATION OF BUOY COORDINATES® INPUT VALUES 
ARE ACTUAL RANGE MEASUREMENTS TAKEN, THE ASSOCIATED A/C 
COORDINATES AT TIME RANGES TAKEN AND AN ESTIMATE OF THE REL- 
ATIVE BUOY ANGLE WITH RESPECT TO A/C HEADING® 

! 

IMPLICIT REAL*8 (A-H,0-Z) 

REAL*8 DSQRT, DCOS, DSI N, DAT AN , P ( 1 0 ) , RP < 1 0 ) , Q ( 1 0 ) , 
1QQ(10> , S (1C) , SS( 10) ,T< 10 I ,R( 10 ) , A (10 ) , B< 10 >, DABS 

1 FORMAT! FlOo A) 

2 F0RMAT(10F8®4) 

3 FORMAT! 1 OP 8o A) 

A FORMAT ( lOFlOo 5 ) 

5 FORMAT ( 'O' T7, * X* ,T17, « Y* ,T24, • ERRORX • ,T34, 'ERRORY* / ) 

6 FORMAT!* *4F10®4) 

INPUT THE ESTIMATE OF THE RELATIVE ANGLE OF THE BUOY® 

READ ( 5 , 1 )THETA 
N = 3 
Z = N 

INPUT THE A/C COORDINATES ( A , B ) AT EACH RANGE MEASUREMENT 0 

READ(5,2>( A! I > , 1 = 1, N) 

READ ( 5 , 3l(B(II,I=l,N) 

READ! 5, A) ( R( I ) , 1=1 , N) 

ABAR=Oo 
BBAR=Oo 
RBAR=Oo 
DC 50 1=1, N 
ABAR=ABAR+A( I ) 

BBAR=BBAR+B ( I ) 

RB AR=RBAR+R ( I ) 

50 CONTINUE 

CALCULATE INITIAL VALUES OF X AND Y (BUOY COORDINATES) USING 
ESTIMATED VALUE OF THE RELATIVE BUOY ANGLE® 

RAD=57 0 2S57T951 
THETAR-THETA/RAD 
X=(RBAR*DCO$(THETAR J + ABAR) VZ 
Y=( RBAR*DS IMTHETAR ) + BBAR )/Z 
WRITE (6,5) 

COMMENCE ITERATIVE PROCESS FOR CONVERGENCE ON VALUES FOR THE 
BUOY COORDINATES® NUMBER OF ITERATIONS MAY BE INPUT HERE® 

DO 100 J=1 , 15 
- V1=0® 

V2=0® 

V3=0® 

V4 = 0® 

V5=0® 

V6=0® 

V7=0o 

DO 70 1=1, N 

P( I ) = ( (X-A! I ) )**2+ (Y-B! I ) )**2>*-71® 5 
PP( I ) = ( ( X-A ( I ) )**2+ (Y-B! I ) )**2 )**o 5 
Q (I ) =( X-A( I ) )**2 
0G( I ) = X- A ( I ) 

S ( I ) = ( Y-B ( I ) )**2 
SS ( I ) = Y- B( I ) 

T ( I ) =( X-A( I ) ) Y (Y-B ( I ) ) 

V1 = V1 + R ( I ) *S ( I ) / P ( I ) 

V2=V2 + R( I ) ^ T ( I ) / P ( I ) 

V3=V3+R( I )*C( I )/P( 1 ) 

V4= V4+R ( I )*G0( I ) /PP(I ) 

V5=V5+R ( I )*SS< I ) /PP ( I ) 

V6=Vfc+ A ( I ) 

V7=V7+B( I ) 

70 CONTINUE 
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D 1 =1 -VI / Z 
D2=V2/Z 
D3 = D 2 
D4=l-V3/Z 

D5=< Y*V2-X*V1+V4+V6 > / Z 

D6=( X*V2-Y*V3+V5+V7 )/Z 

X= ( D5 i -'D4-D3^06 ) / (D1*D4-D3*D2 ) 

Y=(D1*D6-D5*D2)/(D1*D4-D3*D2) 

100 CONTINUE 

WRITE OUT FINAL VALUES OF BUOY COORDINATES AT END OF ITER- 
ATIVE PPOCESSo 

WR I T E ( 6 t 6 ) X » Y 

STOP 

END 
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SONOBUOY LOCATION SENSITIVITY ANALYSIS PROGRAMo DATA OUTPUT 
FROM PROGRAM DETERMINES THF RATE OF CONVERGENCE TO WITHIN 
ACCEPTABLE VALUES OF THE SONOBUOY COORDINATES AS AFFECTED BY 
VARIANCE OF THE FOLLOWING; NUMBER OF RANGE MEASUREMENTS USED 
AND THE AVERAGE RANGE TO THE BUOY, RELATIVE ANGLE OF THE 
SONOBUOY TC THE AIRCRAFT HEADING AND RELATIVE BEARING ESTI- 
MATEo 

IMPLICIT REAL*8 (A-H,0-Z) 

REAL *8 DSQRT,DCOS, DSIN,DATAN,P<10) , PP ( 10 ) , Q ( 10 ) , 
1QQ(10) , S(10 ) , SS( 10) ,T( 10 ) ,R( 101 , A( 10) ,B( 10) 

1 FORMAT ( 3F1 Oo A ) 

2 FORMAT ( 10F8 o 4 ) 

3 FORMAT <10F8 o 4) 

5 FORMAT ( *0*////) 

6 FORM AT ( «0*T5» • THET A • , T 17 , ' X ' , T27, ’ Y 1 ,T37 , »C* ,T47, • D e / ) 

7 FORMAT ( ' 0* 5F1 0 o 4// ) 

8 FORMAT ( • ' 2F10o 4, 1 1 0 ) 

ESTABLISH THE NUMBER OF RANGE MEASUREMENTS TAKENo 

N = 7 
Z = N 

RE AC INPUT VALUES OF A/C COORDINATES FOR EACH OF THE RANGE 
MEASUREMENTSo 

RE AD ( 5 , 2 ) ( A ( I ) , I=1,N) 

READ<5,3)( BID, 1=1, N) 

GO TO 400 
200 WRITE(6,5) 

READ INPUT VALUES OF THE RELATIVE BEARING ESTIMATE AND THE 
ACTUAL BUOY COORDINATES,, 

400 READ(5,1)THETA,C»D 

IFICoEQolOCOJGO TO 300 

ABAR =0 

BBAR=0 

DETERMINE ACTUAL RANGES FROM A/C 10 BUOYo 
DO 50 1=1, N 

R ( I ) =DSQRT ( (C-A( I ) )**2+(D-B( I ) )**2 ) 

AB AR=AB AR+A ( I ) 

B BAR=BBAR + B ( I ) 

50 CONTINUE 

CALCULATE ACCEPTABLE ERROR (TEST) AS o 5% OF RANGE TO BU0Y o 

TESTX=R(2)/200o 
TESTY = TE ST X 
RAD= 57 o 24577951 
RBAR=0 
DO 60 I =1 , N 
RB AR=R B AR + R ( I ) 

60 CONTINUE 

DETERMINE INITIAL VALUES OF BUOY COORDINATES USING ESTIMATED 
VALUE OF RELATIVE BEARINGo 

THET AR=THET A/RAD 
X=(RBAR*DCCS(THETAR )+ABAR) /Z 
Y=(RBAR*DSIN(THETAR)+BBAR) /Z 
WRITE(6,6) 

WRITE(6,7)THETA,X,Y,C,D 
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COMMENCE ITERATIVE PROCESS TO CONVERGE ON ACTUAL VALUES OF 
BUC V COORDINATESo 

DO 100 J =1 i 20 
V1=0 

V 2 = 0 
V3 = 0 
VA=0 
V5=0 

V 6=0 
V7=0 



CALCULATE COEFFICIENT MATRIX AND CONSTANTS FOR SYSTEM OF 
NORMAL EQUATIONS USED IN THE ITERATIVE PROCESSo 

DC 70 1=1, N 

P(I)=< ( X-A( I) )**2+ (Y-B< I ))**2)**lo5 
P P ( I ) = ( (X-A(I) )**2+(Y-B(I ) )**2)**o5 
Q ( I ) =(X-A( I ) ) **2 
QQ( I ) =X-A( I ) 

S ( I ) = ( Y-B C I ) )**2 
S S ( I ) =Y-B ( I ) 

T (I ) = ( X - A ( I ) )*(Y-B( I ) ) 

VI = V 1 +R ( I ) * S ( I ) / P ( I ) 

V2 = V2+R ( I )*T( I ) / P ( I ) 

V3 = V3+R ( I ) *C( I ) / P ( I ) 

V4= V4+R ( I )*GQ( I ) /PP ( I ) 

V5=V5+R ( I )*SS( I ) /PP(I ) 

V 6=V6+A( I ) 

V7=V7+B ( I ) 

70 CONTINUE 
Dl= 1-V1/Z 
D2=V2/Z 
D3=D2 
DA=1-V3/Z 

D5=( Y*V2-X#V1+V4+V6 ) /Z 
D6=( X*V2-Y*V3+V5+V7)/Z 

USE KRAMER'S RULE WITH RESPECT TO THE COEFFICIENTS CALCULAT- 
ED ABOVE TO DETERMINE VALUES OF BUOY COORDINATESo 

X=(D5*DA-D3*D6)/ (D1*D4-D3*D2 ) 
Y=(D1*D6-D5*D2)/(D1*D4-D3*D2) 

WRITE(6,8) X,Y, J 

IF ERROR IS WITHIN ACCEPTABLE VALUE, STOP ITERATIVE PROCESS 
AND NOTE THE NUMBER OF ITERATIONS REQUIRED TO CONVERGE AND 
THE FINAL BUOY COORDINATESo 

TESTXX=DABS(C-X) 

TESTYY=DABS(D-Y) 

IF(TESTXXoGT o TESTX!G 0 TO 100 
I F ( TESTYYo GToTESTYJGO TO 100 
GO TO 200 
100 CONTINUE 
GO TO 200 
300 STOP 
END 
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GROUND SPEED DETERMINATION PROGRAM., FOR USE IN TH F ACTUAL 
DETERMINATION OF GROUND SPEED OF THE AIRCRAFTo INPUTS 
SUP PL I E C ARE TIME BETWEFN RANGE MEASUREMENTS AND THE 
ACTUAL RANGES MEASUREDo 



IMPLICIT RE£L *8 (A-H,0-Z» 

REALMS DSQR T t R ( 1 0 ) 

1 FORMAT ( 10F8o4 ) 

2 FORMAT! 'O' F12o4) 

ESTABLISH NUMBER OF RANGES TAKEN AND TIME B FT WEEN RANGES IN 
SECCNDSo 

N = 3 
Z = N 

T I ME=3o 

INPUT MEASURED RANGESo 

RE AD ( 5 » 1 ) (R(I) ,I=1,N) 

IF(NoGTo3)GC TO 8 

GS = DSQPT ( ( R ( 1 )**2-2o*R(2)**2+R (3)**2 )/2o ) 

GO TO 20 

8 I F ( No GTo A ) GC TO 9 

GS = DSQRT ( ( R ( 1 >** 2-R ( 2 ) **2-R ( 3 ) * *2 + R ( 4 )**2 » /4 0 > 

GO TO 20 

9 I F ( No GTo 5 ) GC TO 10 

G S=DSQRT ( ( 2o*R( 1 »**2-3 0 *R ( 2 > **2 + R ( 3 )**2-R( 41**2+ 

1 R ( 5 ) **2 ) /8o ) 

GO TO 20 

10 I F ( No GT 0 6 ) GC TO 11 

GS=DSQRT ( ( 2o*R(l )**2-2o*R( 2 ) **2-P ( 3 ) **? + R ( A ) **2 
1-R C 5 )**2+R(6)**2)/12> ) 

GO TO 20 

11 I F ( NoGTo 7 ) GC TO 12 

GS=DSQBT(( R(1 ) **2-R ( 2 ) **2+R ( 3 > **2- 2o *R ( 4 ) ** 2+R ( 5 ) ** 2 
1-R(6)**2+R(7)**2 )/12o ) 

GO TO 20 

12 GS- DSQRTU R ( 1 )**2-R ( 2 » ** 2+ R ( 3 I** 2-R ( A )**2-R ( 5 1**2 
1+R(6}**2-R( 7)**2+R (8»**2)/16o ) 

20 GN0SPD=GS*18o/TIME 
WRITE(6,2)GNDSPD 
SI OP 
END 
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GROUND S PR C D DETERMINATION SENSITIVITY ANALYSIS PROGRAM? 
CATA OUTPUT RELATES TO GROUND SPEED ERROR AS AFFECTED BY 
VARIANCE OF RANGE ERRORS, RELATIVE BUOY ANGLE, RANGE TO 
BUOY, NUMBER OF RANGES MEASURED AMD TIME BETWEEN RANGE 
MEASUREMENTSo 

i 

I 

IMPLICIT REAL*8 (A-H,0-Z! 

REAL*8 DSQRT,DABS, R(IO ) ,A( 10), B ( 1 0 ) , ERR ( 10 ) 

1 FORMAT ( 2F10o4l 

2 FORMAT (SFlOoAl 

5 FORMAT ( 'O' T1 , • GROUNDSP EED' T14, ’ERROR' ) 

6 FORMAT ( • • 2F10o 4 ) 

ESTABLISH ACTUAL DISTANCE FLOWN IN NAUTICAL MILES BETWEEN 
RANGE MEASUREMENTS, TIME BETWEEN RANGE MEASUREMENTS AND 
ACTUAL GROUND SPEED TO BE USED AS A TEST FOR DATA GROUND- 
SPEEDo 

E = o 2 5 
T I ME=4o 5 

TEST=E*3600o /TIME 
N = 6 
Z = N 

INPUT ACTUAL BUOY COORDINATES TO SIMULATE VARYING BUOY 
ANGLES AND RANGES RELATIVE TO THE AIRCRAFTo 

RE AD ( 5 , 1 )C , D 

INPUT AIRCRAFT COORDINATES AND RANGE ERROR APPLIED TO 
ACTUAL RANGE TO SIMULATE MEASURED RANGE WITH ERRORo 

RE AD ( 5 , 2 ) < A ( I » , I =1 , N > 

RE AD ( 5 , 2 ) ( B ( I > ,I=1,N) 

WRITE! 6, 5) 

DO ICO J=1 , 70 

READ ( 5 , 2 ) (ERR( I) , I =1 , N ) 

CALCULATE MEASURED RANGE" AS DIFFERING FROM ACTUAL RANGE 
BY ERROR APPLIEDo 
DO 200 1=1 ,N 

R ( I ) =DSQRT ( (C-A(I) )**2+(D-B(I) )**2) 

R ( I )=R( I ) +ERR ( I ) 

200 CONTINUE 

CALCULATE VALUE OF GROUND SPEED FROM INPUT DAT Ao 
I F ( NoGTo 2 ) GC TO 30 

_ GS = DSQRT ( ( R ( 1 )**2-2o*R ( 2 )**2+R ( 3 )**2 »/2o ) 

GO TO 50 

30 I F ( NoGTo 4 I GO TO 31 

G$ = DSQRT( ( R(1 )**2-R(2)**2-R< 3 ) 2+R ( 4 ) /4o ) 

GO TO 50 

31 I F ( No GT 0 5 ) GC TO 32 

GS=DSQRT ( < 2o*R< 1 )**2-3o*R( 2)**2 + R( 3 )**2-R ( 41**2 
1 + R ( 5 )**2 )/8o ) 

GO TO 50 

32 I F ( NoGTo 6 ) GG TO 33 

G S = D SORT ( ( 2o*R(l)*^2-2o*R(2)**2-R(3)**2 + R(4)**2 
l-R(5)**2+R(6)**2)/12o ) 

GO TO 50 

33 I F ( NoGT o 7 ) GO TO 34 

GS = DSQRT ( ( R(1 I** 2-R ( 2 ) **2+R ( 3 I **2- 2- *R < 4 ) **2 + R ( 5 ) **2 
l-R(6)**2+R<7)**2)/12o ) 

GO TO 50 

34 GS = DSORT ( ( R ( 1 ) 2-R ( 2 I 2+ R ( 3 ) 2- R ( 4 ) ** ? - P ( 5 ) "^*2 

1 + R(6)**2-R(7)**2+R(8)*T2 >/16o ) 

50 GNDSPD=GS^36C0o /TI ME 



TEST CALCULATED GRCUND SPEED AGAINST ACTUAL GROUND SPFED TO 
DETERMINE ERRORo 



ERROR=DABS( TEST-GNDSPD ) 
100 WRITE(6,6»GNDSPD, ERROR 
STOP 
END 
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